1 Abstract

The increasing size of the wild boar population causes ecological impacts in natural environments and tremendous economic losses in agriculture. The research objective is to understand if there are changes in the movement patterns during the hunting season. Such knowledge could potentially be used to reduce ecological impacts or improve hunting strategies. Associations between hunting season and the spatial as well as temporal distribution are used together with activity patterns such as traveled distances and speed patterns. Moreover, to detect hunting season associations the overall seasonal patterns are investigated. With density heatmaps the spatial and temporal distribution is investigated. For the speed patterns and traveled distances investigation the trajectory segmentation is applied. Speed and travel distance are investigated by month and day phases with the associated hunting pressure. The key results are that the association between hunting and wild boar activity remains difficult to assess. Nevertheless diurnal patterns are detected that are also found in the literature: wild boar activity is higher during the night than throughout the day.

Keywords: GPS, wild boar, hunting, density heatmaps, animal speed

2 Introduction

Wild Boar (Sus scrofa L.) populations have steadily increased in size across Europe since the 1980s (Massei 2015). Reasons include a high reproductive rate of up to 2.0 (Bieber and Ruf 2005; Keuling, Baubet, and Duscher 2013), limited natural mortality, which can be caused by the wolf (Canis lupus) (Nores, Llaneza, and Alvarez 2008), diseases (Rossi 2011), or starvation due to extreme weather conditions (Okarma 1995). Global warming is expected to further intensify this trend (Melis et al. 2006). This leaves hunting, and to a lesser extent road traffic accidents, as the main contributors to wild boar mortality (Keuling, Baubet, and Duscher 2013; Nores, Llaneza, and Alvarez 2008; Toiga 2008).

The wild boar population did not only expanded in density but also in geographical range (Segura, Acevedo, and Rodriguez 2014). Migration into new areas also occurred due to populations which have been mistakenly reinstated by individuals escaping from farms (Macdonald 2001). Besides causing ecological impacts in natural environments, expansion into territories where animals were not previously found can cause tremendous economic losses in agriculture and forestry (Pimentel et al. 2001). Fischer, Gourdin, and Obermann (2004) stress that the damage increased dramatically since the 1990s. The wild boar, the fifth largest ungulate species in Europe, is an animal that has a far-reaching influence on both natural and anthropogenic areas (Niethammer and Krapp 1986).

The impact of hunting on the wild boar is often indirectly tracked by patterns of wild boar related damage to the flora. et al. Thurfjell H. (2009) added GPS data to his analysis, which is, however, not focussing on the influence of the hunting season and is also limited to Sweden. Since seasonal variations are expected to be influenced by temperature, the findings of this study cannot easily applied to Switzerland. A study from Lower Saxony in Germany found hiding and running away as the two main responses to disturbance by hunting of the wild boars (Sodeikat and Pohlmeyer 2003). However, drive hunts were investigated in this study which has only limited application in Switzerland, where hunting is dominated by blind hunting.

As a result there is currently limited knowledge on the influence of hunting on movement patterns of wild boar. Such knowledge could potentially be used to reduce ecological impacts or improve hunting strategies.

3 Research Questions

This paper intense to address the lack of understanding of how hunting affects the behaviour and the movement patterns of wild boar. Specifically the focus is to investigate the following research questions (RQ) and hypotheses (H).

  • RQ 1: Is there an overall shift in active areas when hunting season starts or during hunting season?

    • H 1.1: The wild boars activity area (spatial and temporal distribution) is different during and before/after hunting season.
    • H 1.2: During the hunting season the diurnal (spatial, temporal and activity) rhythms differ during and before/after hunting season diurnal rhythms.
  • RQ 2: Are there short or medium term effects related to hunting similar to scare off measures?

    • H 2.1: Hunting leads to higher speeds and higher covered distances due to animal fleeing.
    • H 2.2: The activity levels are higher due to (constant) hunting disturbances.

The implemented analysis and code can be displayed by clicking the code button on the right side of the analysis. First, the available data is explored and investigated. Next step includes the description of the used methods in this paper. Afterwards, the results are presented and are further discussed in the last section. Also, the limitations are listed.

4 Materials and Methods

4.1 Study Area

The study area is located at the north east bank of Lake Neuchâtel (N 1197562 - 1207630, E 2564229-2577023 CH1903+ / LV95 (EPSG 2056)), which touches three Swiss cantons: Bern, Neuchâtel and Vaud. The area is covered by four main habitat types: agricultural land, forest, open areas like pastures or meadows, and wetlands. The study area contains protective areas.

Specific hunting activities are not recorded in this area since it is regulated by seasonal permits. Hunting season typically lasts from August until January with each canton having specific regulations regarding the details.

4.2 Position data

The used data set contains the positions of 19 different wildboar, tagged between 2014-05-28 21:01:14 UTC and 2016-10-18 02:00:43 UTC in the context of a research project at ZHAW investigating different repelling systems to protect land crop from wild boar (ZHAW, n.d.). GPS collars were used which typically recorded the position of the animals every 15min. In some cases the collars had to be replaced during the study period. To nonetheless be able to properly assign the data to the correct animal, all wild boar received names. Further details can be found in Section Data.

Additionally, the research project also provided information regarding the different habitats and crops in the study area. This data set is limited in the sense that it has no temporal component as well. As a result, each spatial area is tagged with one type of habitat or crop, although especially agricultural areas might change throughout the year.

4.3 Methods

Different methods were explored to answer the RQs: seasonal density heatmaps to investigate the seasonal impact of the hunting season on the spatial distribution, calculate traveled distances by animal over time to examine the changes of habitats and analysis of speed distribution.

Density heatmaps are used to answer RQ 1. Within spatial statistics, the visualization of data with their basic landmarks and geographic context is essential (Kahle and Wickham 2013). With the density heatmaps seasonal differences can be investigated. Density heatmaps are made for each month and day phase.

In order to answer the RQ 2 and 3 it is necessary to have the information about the different trip segments. (Laube and Purves 2011) introduced one segmentation method. With static fixes as “those whose average Euclidean distance to other fixes inside a temporal window v is less than some threshold d,” the trajectories segmenting is done. However, this method is applied and considered that the data set consists of more than one animal. The applied analysis allows the interpretation and comparison between different months and seasonal characteristics.

5 Data

5.1 Dataset

The used data set contains the positions of 19 different wildboar, tagged between 2014-05-28 and 2016-10-18 in an area at the north east bank of Lake Neuchâtel

Table 1 shows all the available different wild boar animals of the data set. It is possible to investigate the different variables interactively. Moreover, the GPS trajectories of all animals are available. In Table 2 a simple subset of the data is shown. A consistent sampling rate is not available.

Table 1: Wild boar metadata interactive searchable.

DT::datatable(wildschwein_metadata)
  • TierID: The unique identifier for the animal
  • TierName: Name of the ANmial
  • CollarID: The used collars have their own ID
  • Datetime
  • Coordinates
  • day: day phases
  • moonilumination: Important to define full moon occurances.

Table 2: Available GPS trajectories and the given data structure.

# Include date information
wildschwein_BE <- wildschwein_BE %>%
  mutate(month=month(wildschwein_BE$DatetimeUTC)) %>%
  mutate(daymonth=day(wildschwein_BE$DatetimeUTC))

wildschwein_BE$month <- month.name[wildschwein_BE$month]

head(wildschwein_BE)
## # A tibble: 6 x 10
##   TierID TierName CollarID DatetimeUTC                      E       N day       
##    <int> <chr>       <int> <dttm>                       <dbl>   <dbl> <chr>     
## 1      1 Ueli        12272 2014-05-28 21:01:14.000000  2.57e6  1.20e6 Tag       
## 2      1 Ueli        12272 2014-05-28 21:15:18.000000  2.57e6  1.20e6 Abenddaem…
## 3      1 Ueli        12272 2014-05-28 21:30:13.000000  2.57e6  1.20e6 Abenddaem…
## 4      1 Ueli        12272 2014-05-28 21:45:11.000000  2.57e6  1.20e6 Abenddaem…
## 5      1 Ueli        12272 2014-05-28 22:00:33.000000  2.57e6  1.20e6 1Nachtvie…
## 6      1 Ueli        12272 2014-05-28 22:15:16.000000  2.57e6  1.20e6 1Nachtvie…
## # … with 3 more variables: moonilumination <dbl>, month <chr>, daymonth <int>

5.2 Data Exploration

For various reasons, not all animals were tracked over the entire time. In Figure 1 it is investigated for which months data is available per animal. This is cross-correlated with the hunting season to make sure that only animals with data from both periods, hunting and no hunting, is included in the study.

par(mfrow=c(1,2)) 

wildschwein_BE %>%
  group_by(TierID) %>%
  ggplot(aes(x=DatetimeUTC, y=TierName)) +
  geom_line() +
  theme_bw() + ylab("Wild Boar Name") + xlab("Timeline")


wildschwein_BE %>%
  group_by(TierID) %>%  
  ggplot(aes(x=month(wildschwein_BE$DatetimeUTC), y=TierName)) +
  geom_point() +
  scale_x_continuous(breaks=seq(1,12,1))+
  theme_bw() + ylab("Wild Boar Name") + xlab("Months")
Figure 1: Left) Overview of all available GPS trajectories for each wild boar. Right) Overview of the available month by each boar.Figure 1: Left) Overview of all available GPS trajectories for each wild boar. Right) Overview of the available month by each boar.

Figure 1: Left) Overview of all available GPS trajectories for each wild boar. Right) Overview of the available month by each boar.

5.3 Defining Subsample

Hunting season is roughly from August to January, with the main season from November to January (details see Section 6.1. As such, we will include only animals where we have date covering the main hunting season as well as at least two months outside of the hunting season. As a result, we will use the data from:

  • Sabine, Miriam, and Caroline: covering January - December (12 months)

  • Amos: covering July to April (10 months)

  • Ruth: covering November until July (9 months)

  • Rosa: covering November until June (8 months)

Except for Amos, all of these animals are female. As a result, we will not have sufficient statistics to study potential differences in the behaviour of male and female wild boar.

5.4 Spatial and Temporal Distribution

With the investigation of the temporal data availability and the discussed conditions in Section 4.3 the subsample can be defined. Accordingly, the spatial and temporal distribution can be investigated to get an overview of the available data. Moreover, it is possible to interact with the leaflet application and to compare different months. It is observable that the wild boars move between Lake Neuchâtel and Lake Murten.

# Generate subsample
ws_final <- wildschwein_BE [wildschwein_BE$TierName %in% c("Sabine", "Miriam", "Caroline", "Amos", "Ruth", "Rosa"), ] 

# Reproject data to use it for leaflet
ws_final_lv95 = st_as_sf(ws_final, coords = c("E", "N"), crs = 2056)
ws_final_wgs84 = st_transform(ws_final_lv95, 4326)

# We create an SF object to cast points to linestring
ws_final_wgs84 <- ws_final_wgs84 %>% 
  dplyr::group_by(TierName, month) %>% 
  dplyr::summarise() %>%
  st_cast("LINESTRING")

# The SF object is splitted to a list by month 
wildboar.df <- split(ws_final_wgs84, ws_final_wgs84$month)

# Generate color palette for leaflet
factpal <- colorFactor(topo.colors(6), ws_final_wgs84$TierName)

# Leaflet map with baseGroups to get an overview
boarMap <- leaflet() %>% addTiles()
names(wildboar.df) %>%
  purrr::walk( function(df) {
    boarMap <<- boarMap %>%
      addPolylines(data = wildboar.df[[df]],
                   label=~as.character(TierName),
                   popup=~as.character(TierName),
                   group = df,
                   color = ~factpal(TierName),
                   fillOpacity = 1,
                   smoothFactor = 0.2
                   )
  })

boarMap %>%
  addLayersControl(
    baseGroups = names(wildboar.df),
    options = layersControlOptions(collapsed = FALSE)
  )

6 Results

6.1 Defining Hunting Season

The habitat of the wild boars lies within three different political districts and each district has its own hunting regulations. In Table 4 a meta-analysis of the different hunting regulations is shown. Hunting season will be marked as none, low, medium, or high, depending on month and time of day. The following table shows the exact details of these regulations together with details how this was modeled in this analysis. Noteworthy is that full moon is a key element of defining the hunting season and is further investigated in Section 6.2.

Table 4: Hunting regulations in relevant cantons including indicator for hunting season (HS) and which day times this indicator will receive in the data.
Month Bern Vaud Fribourg Indicator HS time of day
Feb - July - potentially stalking - none all
Aug outside forest day only: Mon / Sat stalking day only: Mon, Tue, Thu, Fri, Sat - low day
Sep Mon-Sat day only all day only: Mon, Tue, Thu, Fri. Stalking: +Sat Mon, Tue, Thu, Sat day only. First half only outside forest. No protected areas medium day
Oct Mon, Wed, Sat day only - Mon-Sat, protected areas: 2nd half Mon, Tue, Thu, Sat. Day only low day
Nov

Mon, Wed, Sat day only,

2nd half: 10 days around full moon

all day only: Mon, Tue, Thu, Fri, Sat all areas day only: Mon-Thu, Sat high

day

around full moon: all

Dec Mon, Wed, Sat day only, 10 days around full moon all day only: Mon, Tue, Thu, Fri, Sat all areas day only: Mon-Thu, Sat high

day

around full moon: all

Jan Mon, Wed, Sat day only, 10 days around full moon all day only: Mon, Tue, Thu, Fri, Sat day only: Mon-Thu, Sat
some limitations in protected areas
high

day

around full moon: all

With the gathered information about hunting regulations the GPS trajectories are enriched with the information about the hunting exposure.

ws_final <- ws_final %>%
  mutate(hunting = case_when( (month=="August" | month=="October") & day=="Tag"  ~  "low",
                              month =="September" & day=="Tag"  ~  "medium",
                              (month=="November" | month=="December" | month=="January") & 
                                (day=="Tag" | moonilumination>0.75) ~ "high",
                              TRUE ~ "none"))

ws_final <- ws_final %>% mutate(month = factor(month,
                    levels = c("January", "February", "March", "April",
                               "May", "June", "July", "August",
                               "September", "October", "November", "December"),
                    ordered = TRUE))

ws_final <- ws_final %>% mutate(day = factor(day,
                    levels = c("Abenddaemmerung","1Nachtviertel","2Nachtviertel",
                               "3Nachtviertel","4Nachtviertel","Nacht","Morgendaemmerung","Tag"),
                    ordered = TRUE))

6.2 Full moon detection

From November to January hunting can be done in the nights 6 days before and 4 days after full moon. The data includes a information on moon illumination that can be used to determine when full moon occurs and select the relevant days.

A simple cut as used above on “moon illumination” is giving exactly the required 10 days around full moon. However, filtering the correct days in a somewhat automated way would require a much more sophisticated algorithm (e.g. peak finder on distribution of “moon illumination” to determine full moon or pairing the data with a calendar). The gain in precision was found to be smaller than the overall uncertainties and variations in this study and therefore decided that the simplified cut is sufficient.

Table 3: Full moon day for each month detected with moon illumination >0.75.

fullmoondays <- ws_final %>%
  filter(moonilumination>0.75) %>%
  group_by(month) %>%
  summarise(fullmoonday=n_distinct(daymonth))

DT::datatable(fullmoondays)

6.3 Visited Habitats and Crops

The hunting season might influence the preferred habitat or crop of the wild boar. To investigate this question, the visits of each animal to the different habitats and crops is correlated with the hunting season and its different intensities. This is shown in Figure 2.

Since seasonality in general and availability of the different crops can also have an influence the number of visits per animal and month are also studied. This is shown in Figure 3. Both plots only show the most frequently visited habitats and crops. It has to be pointed out that “others” also include those areas where no habitat or crop information was available (“n/a”).

The first thing that stands out is that each animal has different preferences and that there are no obvious commonalities between all of them. There is also no clear influence of the hunting season since low and high hunting activities tend to have similar patterns. The difference in visited crops and habitats when no hunting takes place can mostly be explained by availability of crops in spring & summer that are not or less available during fall and winter.

fields <- read_sf("./Feldaufnahmen_Fanel.gpkg")

# Join wild boar data set with field information based on coordinates
ws_filter_st <- st_as_sf(ws_final, coords = c("E", "N"), crs = 2056)
ws_join <- st_join(x=ws_filter_st, y=fields)

# Check visited habitats in different hunting seasons
ws_sum_hunt <- ws_join %>%
  st_drop_geometry() %>%
  mutate(Frucht = ifelse(is.na(Frucht),"other",Frucht),
    Frucht = fct_lump(Frucht, 7,other_level = "other"),
    ) %>%
  group_by(TierName, Frucht, hunting) %>%
  summarise(n=n()) %>%
  group_by(TierName, hunting) %>%
  mutate(n_tot=sum(n), 
         n_rel=n/n_tot) %>%
  ungroup() %>%
  mutate(
    Frucht = fct_reorder(Frucht, n,sum, desc = TRUE)
  )

# Polar plot for results
ggplot(ws_sum_hunt, aes(hunting, n_rel, fill=Frucht)) +
  geom_col() +
  coord_polar() +
  facet_wrap(~TierName)  +
  theme_bw() + labs(x="Month", y="Number of visits")
Figure 2: Visited habitats for the different hunting levels and for each animal. Agriculture is detailed by crop type. Shown are only the most important habitats and crops.

Figure 2: Visited habitats for the different hunting levels and for each animal. Agriculture is detailed by crop type. Shown are only the most important habitats and crops.

# Check visited habitats each month
ws_sum_month <- ws_join %>%
  st_drop_geometry() %>%
  mutate(
    Frucht = ifelse(is.na(Frucht),"other",Frucht),
    Frucht = fct_lump(Frucht, 7,other_level = "other"),
    ) %>%
  group_by(TierName, Frucht, month) %>%
  summarise(n=n()) %>%
  group_by(TierName, month) %>%
  mutate(n_tot=sum(n), 
         n_rel=n/n_tot) %>%
  ungroup() %>%
  mutate(
    Frucht = fct_reorder(Frucht, n,sum, desc = TRUE)
  )

# Polar plot for results
ggplot(ws_sum_month, aes(month, n_rel, fill=Frucht)) +
  geom_col() +
  coord_polar() +
  facet_wrap(~TierName) +
  theme_bw() + labs(x="Month", y="Number of visits")
Figure 3: Visited habitats for each month and animal. Agriculture is detailed by crop type. Shown are only the most important habitats and crops.

Figure 3: Visited habitats for each month and animal. Agriculture is detailed by crop type. Shown are only the most important habitats and crops.

6.4 Seasonal and Diurnal Density Heatmaps

After the reprojection and transformation it is possible to apply the ggmap function (ggmap cannot consume LV95). In Figure 2 the density heatmaps for each month including all GPS trajectories of all wild boars are shown. To answer H1.1 the density heatmaps are needed: with the density heatmaps it is possible to get information about the spatial and temporal distribution of the wild boar data for each month. Moreover, plotting each month allows to compare month during and before/after hunting season.

# Reproject data to use it for leaflet
ws_final_wgs84 <- st_transform(ws_final_lv95, 4326)



#'* Extract month for facet wrap *
ws_final_wgs84$day2 <- day(ymd(ws_final_wgs84$DatetimeUTC))
ws_final_wgs84$dayname <-format(ws_final_wgs84$DatetimeUTC, format="%m-%d")
ws_final_wgs84 <- ws_final_wgs84 %>% mutate(month = factor(month,
                    levels = c("January", "February", "March", "April",
                               "May", "June", "July", "August",
                               "September", "October", "November", "December"),
                    ordered = TRUE))

#'* With density maps the seasonal patterns can be detected. *

ws_final_wgs84 <- st_as_sf(ws_final_wgs84, coords = c("lon","lat"))
ws_final_wgs84 <- ws_final_wgs84 %>% tidyr::extract(geometry, c('lon', 'lat'), '\\((.*), (.*)\\)', convert = TRUE)

#'* Define boundingbox (region Neuenburg/Murten and load an stamen basemap*
sbbox <- make_bbox(lon = c(6.973915, 7.129238), lat = c(46.96355, 47.025072), f = .1)
wildRegion <- get_map(location = sbbox, color = "bw", source = "stamen")
wildRegion <- ggmap(wildRegion)

ws_final_wgs84 <- ws_final_wgs84[ -c(9) ]


#'* create density maps for each month*
wildRegion + ggplot2::stat_density2d(aes(x = lon, y = lat, fill = ..level..), data = ws_final_wgs84, alpha = 0.25, bins = 50,  geom = "polygon") +
  scale_fill_gradient(low = "black", high = "red") +
  facet_wrap(~month, ncol = 3) +
  theme_bw() + 
  theme(legend.position="none")
Figure 4: Density heatmaps for each month including all GPS trajectories of all wild boars. Coordinates are transformed into WGS84 due to the fact that the ggmap function cannot consum LV95.

Figure 4: Density heatmaps for each month including all GPS trajectories of all wild boars. Coordinates are transformed into WGS84 due to the fact that the ggmap function cannot consum LV95.

The density heatmaps for each day phase are on key method to answer H1.2: during the hunting season the diurnal (spatial, temporal and activity) rhythms differ before/after hunting season diurnal rhythms. With the density heatmaps the spatial and temporal distribution can be investigated. Plotting by day phase allows to compare the different day stages with each other.

#'* create density maps for day category*
wildRegion + stat_density2d(aes(x = lon, y = lat, fill = ..level..), alpha = 0.25, bins = 50, data = ws_final_wgs84, geom = "polygon") +
  scale_fill_gradient(low = "black", high = "red") + 
  facet_wrap(~ day) +
  theme_bw() + 
  theme(legend.position="none")
Figure 5: Density heatmaps for each day phase including all GPS trajectories of all wild boars.

Figure 5: Density heatmaps for each day phase including all GPS trajectories of all wild boars.

6.5 Traveled Distances and Speed by Animal

In order to answer RQ 2 the traveled distances and speed by animal is investigated. Both parameters might give insights if hunting leads to higher speeds and/or higher covered distances due to animal fleeing. Accordingly, the GPS trajectories are enriched with the speed and traveled distance information. For the analysis the trip segmentation method is applied that is described in the methods section. Moreover, the speed information is calculated out of the distance and the time difference.

# Calculate distances between different steps with step width between 1 - 3 using both, lag and lead. Calculated than the mean of these distances and plot the result as histogram. Save only the distance and also include a marker if the segment is considered static or not.

# split animals into list of lists
ws_final_list = split(ws_final, ws_final$TierName)


rle_id <- function(vec){
  x <- rle(vec)$lengths
  as.factor(rep(seq_along(x), times=x))
}

# Add step mean as attribute
for (i in 1:length(ws_final_list)) {
  
ws_final_list[[i]] <- ws_final_list[[i]] %>%
  mutate(
    stepMean = rowMeans(                       
      cbind(                                   
        sqrt((lag(E,3)-E)^2+(lag(E,3)-E)^2),   
        sqrt((lag(E,2)-E)^2+(lag(E,2)-E)^2),   
        sqrt((lag(E,1)-E)^2+(lag(E,1)-E)^2),   
        sqrt((E-lead(E,1))^2+(E-lead(E,1))^2),  
        sqrt((E-lead(E,2))^2+(E-lead(E,2))^2),
        sqrt((E-lead(E,3))^2+(E-lead(E,3))^2)  
        )                                        
    )
  ) 

  ws_final_list[[i]] <- ws_final_list[[i]]  %>%
  mutate(
    static = stepMean < mean(ws_final_list[[i]]$stepMean,na.rm = TRUE)
  )
  
  # Add segmentation
  ws_final_list[[i]] <- ws_final_list[[i]]  %>%
    mutate(segment_id = rle_id(static))  
    
 # Drop all observation that have a too big time difference
  ws_final_list[[i]] <- ws_final_list[[i]]  %>%
  group_by(segment_id) %>%
  mutate(duration = as.integer(difftime(max(DatetimeUTC),min(DatetimeUTC),"mins"))) %>%
  filter(duration < 20) %>%
  ungroup()
  
}

#revert to single list
ws_final <- do.call(rbind, lapply(ws_final_list, data.frame)) 

# Calculating Speed per Animal
ws_final <- ws_final %>% 
  group_by(TierName) %>%
  mutate(speed = stepMean/duration *3600/1000)   # Speed in km/h

In order to get insights about the traveled distance the mean, sum and median traveled distances are calculated for each day phase. In Table 5 It is possible to sort the variables in ascending or descending order. Due to the lack of information about the definition of the day phase (e.g., their duration) the sum could not be properly normalized. As such this indicator has to be carefully with this in mind. With mean and median traveled distances it is possible to get insights about the overall distribution of the traveled distances by day phase between two observations. The mean and median traveled distances are significantly higher than during the day. Only at dawn the values are even smaller than during the day.

Table 5: Day phases with mean, sum and median traveled distances.

#'* Investigate traveled distance during the day categories *
ws_final_distance <- ws_final %>%
  group_by(day) %>%
  summarise(         
    mean_distance = round(mean(stepMean,na.rm=T),2),
    sum_distance = round(sum(stepMean,na.rm=T),2),
    median_distance = round(median(stepMean,na.rm=T),2)

  )

DT::datatable(ws_final_distance)

Mean traveled distance combined of all wild boars for each month and day phase are shown in Figure 4. The traveled distance gives insights about high and low activity levels. These activities are also found in the literature as H. Thurfjell, Spong, and Ericsson (2014) states higher activities during the night and lower activates throughout the day, in agreement with our findings.

#'* Investigate covered distance during the day categories *
wildDay <- ws_final %>%
  group_by(month, daymonth, day) %>%
  summarise(         
    mean_distance = mean(stepMean,na.rm=T)
  )

ggplot(wildDay, aes(x = daymonth, y = mean_distance)) + 
  geom_line(aes(color = day, linetype = day)) + 
  scale_color_manual(values = c("dodgerblue", "cornflowerblue", "darkorchid", "mediumorchid" , "turquoise", "darkorange", "goldenrod", "brown")) +
  facet_wrap(~ month) +
  theme_bw() +
  labs(x ="Days [d]", y = "Mean traveled distance [m]")
Figure 6: Median walking distance of all tracked wildboars for each month.

Figure 6: Median walking distance of all tracked wildboars for each month.

Activity levels seem to be lower between mid-April until end of June that could potentially be related to less hunting activity or the fact, that food is easier available due to agricultural activities. This is unclear from the available data. Most striking is, however, the dip in activity for about a week in September. The authors have no explanation for this drop in covered distance.

The variability during the day phases is more visible with the boxplot comparison. During the day the activity is low and increases during the different night stages. All blue and violet lines represent activities during the night. Noteworthy is the overall activity low during the mid of September. Highest activities are recorded during November and December.

ws_final %>%
  ggplot( aes(x=day, y=stepMean, fill=day)) +
    geom_boxplot() +
    scale_fill_viridis(discrete = TRUE, alpha=0.6) +
    theme_bw() +
    theme(
      legend.position="none",
      plot.title = element_text(size=11)
    ) +
  labs(x="Day phase", y="Median traveled distance [m]") + theme(axis.text.x = element_text(angle = 45, hjust = 1))
Figure 7: Median walking distance of all wildboars at the different diurnal phases.

Figure 7: Median walking distance of all wildboars at the different diurnal phases.

Figure 8 shows how the covered distances relate not only to diurnal and seasonal phases but also the different hunting levels. The plot shows total values of the median traveled distance per month and hunting season together with the medians per hunting level over all months to highlight differences throughout the year. Again, there are striking differences per animal with Caroline for example showing significantly larger distances traveled during high hunting activities while Sabine and Ruth are quite the opposite. The reason is currently unclear and requires further investigations.

# Calculate median per month and hunting season
ws_dist_month <- ws_final %>%
  filter(!static) %>%
  group_by(hunting, TierName, month) %>%
  summarise(
    mean_dist   = mean(stepMean, na.rm=TRUE),
    sum_dist    = sum(stepMean, na.rm=TRUE),
    median_dist = median(stepMean, na.rm=TRUE)
  )
# Calculating median per hunting season
ws_dist_year <- ws_final %>%
  filter(!static) %>%
  group_by(hunting, TierName) %>%
  summarise(
    mean_dist   = mean(stepMean, na.rm=TRUE),
    median_dist = median(stepMean, na.rm=TRUE)
  )
# Plot results
ws_dist_month %>%
  ggplot(aes(month, median_dist)) +
  geom_point(aes(color=hunting)) +
  geom_hline(data=ws_dist_year, aes(yintercept=median_dist, color=hunting)) +
  theme_bw() + labs(x="Month", y="Median traveled distance [m]") + theme(axis.text.x = element_text(angle = 45, hjust = 1)) +
  facet_wrap(~TierName) 
Figure 8: Median walking distance of individual wildboars for each month and hunting activity.

Figure 8: Median walking distance of individual wildboars for each month and hunting activity.

Figure 9 shows how the median speed relate not only to diurnal and seasonal phases but also the different hunting levels. The plot shows total values of the median speed per month and hunting season together with the medians per hunting level over all months to highlight differences throughout the year.

Again, there are striking differences per animal and these differences cannot even be related to the traveled distances as shown in Figure 8. For example, Ruth showed lower traveled distances during high hunting activites but also runs faster in this phase. Agian, the reason is currently unclear and requires further investigations.

# Calculate median per month and hunting season
ws_speed_month <- ws_final %>%
  filter(!static & speed < 60) %>%
  group_by(hunting, TierName, month) %>%
  summarise(
    mean_speed   = mean(speed, na.rm=TRUE),
    max_speed    = max(speed, na.rm=TRUE),
    median_speed = median(speed, na.rm=TRUE)
  )
# Calculating median per hunting season
ws_speed_year <- ws_final %>%
  filter(!static & speed < 60) %>%
  group_by(hunting, TierName) %>%
  summarise(
    mean_speed   = mean(speed, na.rm=TRUE),
    max_speed    = max(speed, na.rm=TRUE),
    median_speed = median(speed, na.rm=TRUE)
  )
# Plot results: median
ws_speed_month %>%
  ggplot(aes(month, median_speed)) +
  geom_point(aes(color=hunting)) +
  geom_hline(data=ws_speed_year, aes(yintercept=median_speed, color=hunting)) +
  theme_bw() + labs(x="Month", y="Median speed [km/h]") + theme(axis.text.x = element_text(angle = 45, hjust = 1)) +
  facet_wrap(~TierName) 
Figure 9: Median speed of individual wildboars for each month and hunting activity.

Figure 9: Median speed of individual wildboars for each month and hunting activity.

7 Limitations

The analyzed GPS trajectories in this paper have no consistent sampling rate. However, this rate is not perfectly suitable to analyze associations between hunting and moving patterns such as fleeing, speed variations or traveled distance. The information between two observations can only assumed. Moreover, the data coverage is rather low and the sample size with six animals is too small to get significant results.

The investigation of the wild boar hunting season is difficult to assess due to the fact that the habitat of the wild boars lies within three different political districts. In other words, three different hunting regulations have to be considered and therefore it is difficult to define the hunting season for this area. To improve on this point it would be helpful to correlate the data with actual killings due to hunting or potentially even car accidents in the region. Each of these events would ideally be tracked with position and time.

Furthermore H. Thurfjell, Spong, and Ericsson (2014) stress that the movement speed of wild boar is impacted by daylight, season, and weather. During the night and high temperatures higher movement speed is recorded and snow and precipitation lead to lower movement speed. This analysis didn’t include any weather information.

8 Discussion

This paper intended to provide answers on the association between hunting season and spatial and temporal distribution of wild boars in Switzerland. Furthermore a research goal is to investigate the short or medium term effects related to hunting with focus on speed and traveled distance.

The main hunting season is detected and defined between September and January with highest hunting activities expected to happen between Novemer and January. It was necessary to have wild boar data during and before/after the hunting season. For this purpose the available GPS trajectories are investigated. An overview of all available trajectories for each boar and an overview of the available months by each boar is used to define the wild boar selection. After the data exploration 6 out of 19 animals are used for all further analysis in this paper. Except for Amos, all of these animals are female. As a result, it is not possible to have sufficient statistics to study potential differences in the behaviour of male and female wild boar. According to the leaflet application with the spatial and temporal distribution of the wild boar sample the home ranges are always similar and only different during April.

8.1 Findings regarding Research Question 1

The first research question investigated if there isan overall shift in active areas when hunting season starts or during hunting season.

The first hypothesis was investigated with the monthly density heatmaps. The hypothesis states that the wild boar activity (spatial and temporal distribution) is different during hunting season than before or after hunting season. With the defined hunting season (November - January) there is no spatial or temporal association possible. Rather it is a seasonal pattern visible. Early in the year the wild boars home range is always location near the shore. During spring the wild boars move further inland. In the later months of the year no distribution pattern is detected. A potential explanation can be the availability of certain crop types preferred by the wild boar as the visited crops change throughout the year.

For the second hypothesis that during the hunting season the diurnal (spatial, temporal and activity) rhythms differ during and before/after hunting season diurnal rhythms the density heatmaps for each day phase deliver important insights. As mentioned the activity is lower during the day. This can be verified with the density heatmaps. During the day phase the spatial distribution is denser than during the night. Furthermore the activity ranges during all night phases

8.2 Findings regarding Research Question 2

RQ 2: Are there short or medium term effects related to hunting similar to scare off measures?

H 2.1: Hunting leads to higher speeds and higher covered distances due to animal fleeing.
H 2.2: The activity levels are higher due to (constant) hunting disturbances.

The second research question investigated short and medium term effects related to hunting similar to scare off measures. In this study it could not be confirmed that hunting activities indeed lead to higher speeds and larger covered distances due to the animals fleeing. The findings were quite different for each animal and no overarching pattern could be found. As such, H 2.1 could neither be validated nor falsified. Further investigations with more and more detailed data is required in this case.

Similarly, it could not conclusively be confirmed that activity levels are overall higher due to (constant) hunting disturbances. There is an indication of a drop in activity during April to June but this can also be correlated to availability of certain crop types that potentially simplifies foraging. As such, H 2.2 could neither be validated nor falsified. Further investigations with more and more detailed data is required in this case.

8.3 Conclusion

Conclusively it is apparent that there are seasonal patterns in the activity and the spatial and temporal are present. However, it is difficult to associate these patterns directly to the hunting season. One of the reasons is the fact that it is difficult to associate the activity behaviour with actual hunting events. Moreover, the data sampling with 6 animals is rather low. The dirunal pattern found e.g. by H. Thurfjell, Spong, and Ericsson (2014) that wild boar is overall more active during night than day in Sweden was also found for wild boar in Switzerland with higher activity levels during dusk and the first half of the night than during the day.

9 References

Bieber, Claudia, and Thomas Ruf. 2005. “Population Dynamics in Wild Boar Sus Scrofa: Ecology, Elasticity of Growth Rate and Implications for the Management of Pulsed Resource Consumers.” Journal of Applied Ecology 42 (6): 1203–13. https://doi.org/https://doi.org/10.1111/j.1365-2664.2005.01094.x.
Fischer, C., H. Gourdin, and M. Obermann. 2004. Spatial behaviour of the wild boar in Geneva, Switzerland: testing methods and first results.” Galemys 16: 149–55.
Kahle, D. J., and H. Wickham. 2013. ggmap: Spatial Visualization with ggplot2.” The R Journal 5.
Keuling, O., E. Baubet, and A. et al. Duscher. 2013. Mortality rates of wild boar Sus scrofa L. in central Europe.” Eur J Wildl Res 59: 805–14. https://doi.org/https://doi.org/10.1007/s10344-013-0733-8.
Laube, Patrick, and Ross S. Purves. 2011. “How Fast Is a Cow? Cross-Scale Analysis of Movement Data.” Transactions in GIS 15 (3): 401–18. https://doi.org/https://doi.org/10.1111/j.1467-9671.2011.01256.x.
Macdonald, D. 2001. The new encyclopedia of mammals. Oxford University Press.
Massei, Giovanna et al. 2015. “Wild Boar Populations up, Numbers of Hunters down? A Review of Trends and Implications for Europe.” Pest Management Science 71 (4): 492–500. https://doi.org/https://doi.org/10.1002/ps.3965.
Melis, Claudia, Paulina A. Szafra?ska, Bogumi?a J?drzejewska, and Kamil Barto? 2006. “Biogeographical Variation in the Population Density of Wild Boar (Sus Scrofa) in Western Eurasia.” Journal of Biogeography 33 (5): 803–11. http://www.jstor.org/stable/3838491.
Niethammer, J., and F. Krapp. 1986. Handbuch der Säugetiere Europas. Wiesbaden: Aula-Verlag.
Nores, Carlos, Luis Llaneza, and Angel Alvarez. 2008. Wild boar Sus scrofa mortality by hunting and wolf Canis lupus predation: an example in northern Spain.” Wildlife Biology 14 (1): 44–51. https://doi.org/10.2981/0909-6396(2008)14[44:WBSSMB]2.0.CO;2.
Okarma, H. et al. 1995. The roles of predation, snow cover, acorn crop, and man-related factors on ungulate mortality in Bialowieza Primeval Forest, Poland.” Acta Theriologica. http://rcin.org.pl/ibs/Content/12344/BI002_2613_Cz-40-2_Acta-T40-nr17-197-217_o.pdf.
Pimentel, David, S. McNair, J. Janecka, J. Wightman, C. Simmonds, C. OConnell, E. Wong, et al. 2001. “Economic and Environmental Threats of Alien Plant, Animal, and Microbe Invasions.” Agriculture, Ecosystems and Environment 84 (1): 1–20. https://doi.org/https://doi.org/10.1016/S0167-8809(00)00178-X.
Rossi, Sophie et al. 2011. New Insights on the Management of Wildlife Diseases Using Multi-State Recapture Models: The Case of Classical Swine Fever in Wild Boar.” PLOS ONE 6 (9): 1–9. https://doi.org/10.1371/journal.pone.0024257.
Segura, A., P. Acevedo, and O. et al. Rodriguez. 2014. Biotic and abiotic factors modulating wild boar relative abundance in Atlantic Spain.” Eur J Wildl Res 60: 469–76. https://doi.org/https://doi.org/10.1007/s10344-014-0807-2.
Sodeikat, Gunter, and Klaus Pohlmeyer. 2003. Escape movements of family groups of wild boar Sus scrofa influenced by drive hunts in Lower Saxony, Germany.” Wildlife Biology 9 (4): 43–49. https://doi.org/10.2981/wlb.2003.063.
Thurfjell, et al., H. 2009. Habitat use and spatial patterns of wild boar Sus scrofa (L.): agricultural fields and edges.” European Journal of Wildlife Research 55: 517–23.
Thurfjell, H., G. Spong, and G. Ericsson. 2014. Effects of weather, season, and daylight on female wild boar movement.” Acta Theriologica 59: 467–72.
Toiga, Carole et al. 2008. “Disentangling Natural from Hunting Mortality in an Intensively Hunted Wild Boar Population.” The Journal of Wildlife Management 72 (7): 1532–39. https://doi.org/https://doi.org/10.2193/2007-378.
ZHAW. n.d. Evaluation eines akustischen Wildschwein-Vergrämungssystems .” https://www.zhaw.ch/en/research/research-database/project-detailview/projektid/842/.

Sources

Hunting Season